Từ sóng liên tục đến ma trận rời rạc
Dao động của một sợi dây hoặc màng được điều khiển bởi phương trình sóng:
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
Để tìm nghiệm $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$, chúng ta phải giải hệ hệ Sturm-Liouville:
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
Việc rời rạc hóa toán tử dẫn đến các phương trình ma trận như $Aw = -0.04 \frac{\rho}{p} \lambda w$. Đối với ma trận tam giác $4 \times 4$, $p(\lambda)$ là khả thi. Tuy nhiên, khi lưới được làm mịn hơn ($n$ tăng lên), chúng ta gặp phải hai rào cản:
- Hạn chế Abel-Ruffini: Không tồn tại lời giải đại số cho nghiệm của các đa thức khi $n \ge 5$.
- Tính nhạy do làm tròn: Trong các hệ thống chiều cao, một thay đổi ở chữ số thập phân $10^{-10}$ của một phần tử có thể làm dịch chuyển các giá trị riêng lên tới hàng đơn vị (hiện tượng đa thức Wilkinson).
Sự cần thiết về mặt số học và thư viện chuyên nghiệp
Các thư viện số chuyên nghiệp (IMSL, NAG) tránh sử dụng đa thức đặc trưng nguyên thủy. Thay vào đó, họ dùng các thủ tục lặp để xấp xỉ:
- Thư viện IMSL: Sử dụng bình phương bé nhất tuyến tính, đường cong spline bậc ba và biến đổi Fourier nhanh.
- Thư viện NAG: Sử dụng xấp xỉ đa thức bình phương bé nhất và phù hợp theo nghĩa $l_1/l_{\infty}$.
Khi xấp xỉ các giá trị riêng cho hệ $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$, chúng ta dựa vào phương pháp bình phương bé nhất rời rạc và khám phá lặp thay vì tìm nghiệm.